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We propose a new description of inclusive diffraction in deep inelastic scattering (DIS). The 
diffractive structure functions are expressed in the dipole picture and contain heavy-quark contri- 
butions. The dipole scattering amplitude, a saturation model fitted on inclusive DIS data, features 
a saturation scale Q s (x) larger than 1 GeV for x = 10 -5 . The qqg contribution to the diffractive 
final state is modeled in such a way that both the large— Q 2 and small—/? limits are implemented. 
In the regime Xp < 0.01 in which saturation is expected to be relevant, we obtain a parameter-free 
description of the HERA data with x 2 /p°ints= 1.2. 



I. INTRODUCTION 

Deep inelastic scattering (DIS) is a process in which a virtual photon is used as a hard probe to resolve the small 
distances inside a proton and study its partonic constituents: quarks and gluons that obey the laws of perturbative 
QCD. When probing with a fixed photon virtuality Q 2 » ^qcd-> an d increasing the energy of the photon-proton 
collision W, the parton densities seen by the photon inside the proton grow. Eventually, at some energy much bigger 
than the hard scale, corresponding to a small value of the Bjorken variable x ~ Q 2 /W 2 , the gluon density probed 
I ^ ■ becomes so large that non-linear effects like gluon recombination become important. One enters a non- linear yet 
weakly-coupled regime of QCD called the saturation regime. 

The transition to the saturation regime is characterized by the so-called saturation momentum Q s (x) = Qo x~ x / 2 . 
This is an intrinsic scale of the high-energy proton which increases as x decreases. Qq ~ A.qcd, but as the energy 

1 increases, Q s becomes a hard scale, and the transition to saturation occurs when Q s becomes comparable to Q. The 

higher Q 2 is, the smaller x should be to enter the saturation regime. Part of the DIS events are diffractive, meaning 
that the proton remains intact after the collision and there is a rapidity gap between that proton and the rest of the 
O final-state particles. Such events are expected to be more sensitive to the saturation regime than the inclusive ones. 

Although the saturation regime is only reached when Q s ~ Q, observables are sensitive to the saturation scale 
i-C . already during the approach to saturation [|J when Aqcd^Qs^Q- For inclusive events in deep inelastic scattering, 
this feature manifests itself via the so-called geometric scaling property: instead of being a function of Q 2 /Qq and 
H ■ x separately, the total cross-section is only a function of r — Q 2 /Q 2 (x), up to large values of r. Experimental 
' measurements of inclusive DIS are compatible with that prediction Q. Recently, it was shown [3] that diffractive 
O^l , observables also feature the geometric scaling behaviors expected when approaching saturation. 

In the saturation regime of QCD, contributions to the cross-sections growing like Q s /Q are important. The 
leading-twist approximation of perturbative QCD, in which Q 2 is taken as the biggest scale, cannot account for 
. ■ such contributions, and therefore is not appropriate to describe the small— a; limit of deep inelastic scattering. As 
leading-twist gluon distributions cannot be used to compute cross-sections, the dipole picture of DIS [j| has been 
developed to describe the high-energy limit. It expresses the hadronic scattering of the virtual photon through its 
fluctuation into a color singlet qq pair (or dipole) of a transverse size r~ 1/Q. The dipole is then the hard probe that 
resolves the small distances inside the proton. 

The dipole picture naturally incorporates the description of both inclusive and diffractive events into a common 
theoretical framework @, 0| , as the same dipole scattering amplitudes enter in the formulation of the inclusive and 
diffractive cross-sections. Different saturation parametrizations of the dipole-proton cross-section have been successful 
in describing inclusive and diffractive HERA data. The pioneering work of [8|, [9j triggered several improvements: the 
"DGLAP-improved" model of [1C| allows to include even high— Q 2 data in the fit and the "BK-inspired" model of 
[ll| incorporates features from the QCD non-linear equations. 

In diffractive DIS, when the invariant mass Mx of the diffractive final-state is much smaller than Q, the dominant 
contribution to the final state comes from the qq component of the photon wavefunction. By contrast, if j3~Q 2 /M x <C 
1, then the dominant contributions come from the qqg component, or from higher Fock states, i.e. from the photon 
dissociation. The main goal of this work is to improve the description of the qqg contribution with respect to previous 
analysis: it will be modeled in such a way that both the large— Q 2 and small—/? limits are implemented. 
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FIG. 1: Representation of 7* — p deep inelastic scattering; inclusive (left) and diffractive (right) events are pictured with the 
relevant kinematic variables: the photon virtuality Q 2 , the energy squared of the 7*— p collision W 2 , and in the case of diffraction 
the momentum transfer t and the invariant mass of the diffractive final state M|. 



Including the contributions of heavy quarks in the models has also been a recent concern, as several aproaches 
observed a decrease of the saturation scale to Q s ~ A-qcd, when trying to include the charm quark in their analysis 
[T3.[l3j. This was problematic, however it was recently shown [3] that it is possible to accomodate the model of 11 1 
with heavy-quark contributions and a saturation scale that stays above 1 GeV for x = 10~ 5 , rather than dropping 
to about 500 MeV as is the case in other studies. Our second goal in this paper is to check whether the dipole 
cross-section of this heavy-quark improved saturation model also describes the inclusive diffraction data from HERA. 

The plan of the paper is as follows. In Section II, we recall the QCD dipole picture for inclusive and diffractive DIS in 
terms of the dipole-proton scattering. In Section III, we discuss in more details the case of diffraction and present the 
different components of the model, highlighting in each case the improvements with respect to previous approaches, 
in particular concerning the inclusion of heavy-quark contributions, and the treatment of impact parameter. Section 
IV discusses how to implement the qqg contribution to the diffrative final state to obtain a unified description that 
features both the large— Q 2 and small—/? limits. In Section V, the results of the comparison with the available HERA 
data are presented, and Section VI is devoted to conclusions. 



II. THE QCD DIPOLE PICTURE OF DEEP INELASTIC SCATTERING 

We focus on diffractive DIS: 7*p^ Xp (see Fig.l). With a momentum transfer t<0, the proton gets out of the 
j*—p collision intact, and there is a rapidity gap between that proton and the final state X whose invariant mass we 
denote Mx- We recall that the photon virtuality is denoted Q 2 , and the 7*— p total energy W. It is convenient to 
introduce the following variables: 

/3 = m ~ , Jo > £P = x/13 . (1) 



■ Q 2 + W 2 ' f Q 2 + M\ 

The 7*— p total cross-section 0'7 O / > ~* X * s usually expressed as a function of x and Q 2 , while the diffractive cross-section 
da^jJ* Xp / dj3dt is expressed as a function of /?, xp, Q 2 , and t. The size of the rapidity gap in the final state is ln(l/xp). 



A. The 7* — > qq wavefunctions 



To compute those cross-sections in the high-energy limit, it is convenient to view the process in a particular frame 
called the dipole frame. In this frame, the virtual photon undergoes the hadronic interaction via a fluctuation into 
a colorless qq pair, called dipole, which then interacts with the target proton. The wavefunctions ^{'"^(z, r; Q 2 ) 
describing the splitting of a virtual photon with polarization A into a dipole are well known. The indices a and (3 
denote the spins of the quark and the antiquark composing the dipole of flavor /. The wavefunctions depend on Q 2 , 
the fraction z of longitudinal momentum (with respect to the 7*— p collision axis) carried by the quark, and the 
two-dimensional vector r whose modulus is the transverse size of the dipole. 
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FIG. 2: The QCD dipole picture of deep inelastic scattering. The left diagram represents 7* — p elastic scattering and (via 
the optical theorem) corresponds to formula ([7]). The right diagram represents diffractive scattering (without possible final 
states containing gluons) and corresponds to formula JB). In this case, the final state (indicated by the vertical dashed line) 
is characterized by t — — A 2 and M| = (n 2 + m*)/(z(l — 2)), with A = q + q' and k — (1 — z)q — zq 1 in terms of the quark 
and antiquark momenta q and q'. Via Fourier transformations, q and q' impose different sizes and impact parameters for the 
dipole in the amplitude and the dipole in the complex conjugate amplitude. 



Formulae giving the functions V'a'"^ can ^ e found in the literature (see for instance [HI). In what follows, we will 
need the functions <E>{ which describe the overlap between two wavefunctions for splitting into dipoles of different 
transverse size r and r' : 



(*, r, r'; Q 2 ) = N c £ [4' af3 (z, r'; Q 2 )] * ^{z, r; Q 2 ) 

a/3 



For a transversely (T) or longitudinally (L) polarized photon, these functions are given by 
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In the above, e/ and m/ denote the charge and mass of the quark with flavor / and 

E 2 } = z(l-z)Q 2 + m) . 



(4) 



(5) 



B. The total cross-section a7„f^ x 



Via the optical theorem, the 7*— p total cross-section is related to the elastic scattering of the virtual photon off 
the proton. In the dipole frame, this happens as follows: at a given impact parameter b, the photon splits into a 
dipole with a given size r which scatters elastically off the proton and recombines back into the photon. Therefore 
the overlap function $a which enters in the computation of the total cross-section is 



$ A (z,|r|;Q 2 ) =]T0{(z,r,r;Q 2 ) . 
/ 

For a virtual photon with polarization A, the total cross-section is then given by (see Fig. 2a): 
4^ X {x,Q i )=2 f d 2 r [ dz$p*(z,\v\;Q 2 ) [ d 2 bT q ,(r,b;x) 



(6) 



(7) 



where the function T q q(r,h;x) is the elastic scattering amplitude of the dipole of size r off the proton at impact 
parameter b. It contains the x dependence, reflecting the fact that in our frame, the proton carries all the energy and 
is therefore evolved up to the rapidity ln(l/a:). In the high-energy limit i<1 we are considering here, T qq does not 
depend on z. 



C. The diffractive cross- section daJJf Xp /df3dt 

The diffractive scattering happens as follows. In the amplitude, the photon splits into a dipole of size r which 
scatters off the proton at a given impact parameter b and dissociates into a final state of invariant mass Mx- The 
same happens in the complex conjugate amplitude, except that the dipole size r' and the impact parameter b' are 
different from r and b. Indeed, the final state is characterized by particular values of Mx (or equivalently 0) and t, 
corresponding to particular momenta of the quark and antiquark in the final state. In coordinate space, this imposes 
two different dipole sizes and impact parameters in the amplitude and the complex conjugate amplitude, therefore 
the functions 4> x (z,r,r';Q 2 ) (see (fSJ) enter in the computation of the diffractive cross-section. For a virtual photon 
with polarization A, the diffractive cross-section is given by (see Fig. 2b): 

o 2 * = «lv ffL OIL f\,,n-,^ P ^- r) 



d/3 dt 



(z, r, r'; Q 2 ) / d 2 b d 2 b' e^ h '~ h % q (r, b; x r )T qq (r\ b'; a*) . (8) 



In the above, the differences between r and r' on one hand, and b and b' on the other hand, are related via Fourier 
transformation to 

K 2 f = z(l-z)Q 2 (l-/3)/0-m 2 f , and A 2 = -t . (9) 

Note that now, the proton is only evolved up to the rapidity Ln(l/xp). This is because some of the energy (M x ) is 
carried by the dipole in order to form the diffractive final state. The dipole is evolved up to a rapidity rn(l//3) and 
the proton up to the rapidity ln(/3/x) =ln(l/xp). The relevant high-energy limit in this case is ip<l. 

Note that to write formula |(HJ), we have neglected possible final states containing gluons. This is justified because 
these are suppressed by extra powers of a s . However, if (3 becomes too small, or if Q 2 becomes too large, the dipole 
will emit soft or collinear gluons whose emissions are accomponied by large logarithms ln(l//3) or ln(Q 2 ) which will 
compensate the factors of a s . This will be discussed in more details in Section IV, when we explain how to implement 
the qqg contribution to the diffractive final state, in order to correctly describe both the small—/? and large— Q 2 limits. 



III. THE SATURATION MODEL FOR THE DIPOLE AMPLITUDE T q? 

Using the dipole picture of deep inelastic scattering, we have expressed the total © and diffractive ([5]) cross-sections 
in the high-energy limit in terms of a single object: the dipole scattering amplitude off the proton T qq (r, b;x). It is 
mainly a non-perturbative quantity, but its evolution towards small values of x (or high energy) is computable from 
perturbative QCD. Evolution equations have been established in the leading ln(l/x) approximation [UGH and, at 
least for central impact parameters, one has learned a lot about the growth of the dipole amplitude and the transition 
from the leading- twist regime T qq ^l towards and into the saturation regime T qq <l. 

Let us recall that this transition is characterized by the saturation scale Q s (x), which increases as x decreases. In the 
following, we shall work in the context of the Balitsky-Kovchegov (BK) evolution to describe the r dependence. 
Indeed, this provides a natural explanation for the geometric scaling properties of the data 0, lij. The impact 
parameter dependence of T qq is still an open problem, it cannot be extracted from perturbative QCD and it is usually 
modeled. In formulae, one writes 

T qq (r, b; x) = 5(b) N(\r\Q s (x),x) (10) 

where we have introduced the factorized impact-parameter profile 5(b). In the following, we detail the different 
components of our model: 5(b) and iV(|r|Q s , x). 



A. The impact-parameter profile 5(b) 



When performing the b integration in formula l[7jl . this contributes only to the normalization via a constant 
factor 2 J d 2 b 5(b) = ctq, (of order 25 mb) caracterising the transverse area of the proton. However, in the case of 
the diffractive cross-section (|8|). the b integration gives the momentum transfer dependence. Experimentally, the 
diffractive cross-section decreases exponentially with \t\ as e Bot , where Bjj is the diffractive slope (of order 6 GeV -2 ). 
This is consistent with the Gaussian profile 5(b) — e^ h /( 2Bd \ which then implies uq — AttBd. 



In the literature, the quantities oo and Bp are usually considered unrelated, however as we have shown, a consistent 
treatment of the impact parameter dependence within the dipole picture implies that this is not the case. To 
summarize, one has: 

% ~ e B ^ =► S(b) = e - b2 /(2S D ) ^ 47rBz3 . (n) 

dt 



B. The heavy-quark improved IIM saturation model for N(\r\Q s ,x) 

The Iancu-Itakura-Munier (IIM) saturation model is inspired by universal properties [19[ of solutions of the BK 
equation [l6j |. The most important feature is probably the geometric scaling regime: at small values of x, instead of 
being a function of a priori the two variables r = |r| and x, N is actually a function of the single variable rQ s {x) up 
to inverse dipole sizes significantly larger than the saturation scale Q s {x). If rQ s > 1 then N = 1 and the scaling is 
obvious. We insist that the scaling property is a non-trivial prediction for rQ s -C 1, when N is still much smaller than 
1. 

Of course the geometric scaling window has a limited extension: at very small dipole sizes, deep into the leading- 
twist regime, the scaling breaks down. Universal scaling violations [l9| due to x not being small enough have also 
been derived and are implemented in the IIM model, which is therefore a function of rQ s and x. Recently, a new type 
of geometric scaling violations has been predicted, due to the inclusion of Pomeron loops in the evolution (the 
BK equation only resums fan diagrams) . These violations transform the geometric scaling regime into an intermediate 
energy regime, as they arise at very small values of x in the so-called diffusive scaling regime. This new regime is 
likely out of the reach of HERA and we shall not address it in this study. 

In the IIM model, the saturation scale is parametrized by 



and the dipole amplitude is given by 



N(rQ s ,x) 



Q s ( x ) = (^y GeV 



( I §^) 27C exp 

I _ e -4aln 2 (/3rQ s ) 



21n 2 (rQ s /2) 
k\ \n(l/x) 



(12) 



for rQ s < 2 
for rQ s > 2 



(13) 



with a and (3 uniquely determined from the conditions that N and its derivative are continuous at rQ s = 2. The 
amplitude at the matching point is chosen to be No — 0.7. 

In this work, we shall consider the IIM saturation model [ll| extended in [T3| to include heavy quarks (with 
m c — lA GeV, TOf, = 4.5 GeV, and rrif = 0.14 GeV for the light flavors). The coefficient k = 9.9 is obtained from the 
BFKL kernel while the critical exponent j c — 0.7376 is fitted to the HERA measurements of the proton structure 
function, along with the remaining parameters. The saturation scale parameters are A = 0.2197 and xq = 1.632 10~ 5 
and the cross-section at saturation is (Tq = 70.26 GeV -2 (or 27.36 mb). Note that, via UQ — AirBr), this corresponds to 
the diffractive slope Bjj = 5.591 GeV" 2 , which is in agreement with the experimental observations [2ll I22I]. 



C. The qq components of the diffractive structure functions 

Let us introduce the transverse and longitudinal diffractive structure functions F®' 3 (f3, xp, Q 2 ) and F®' 3 ((3, xp, Q 2 ). 

They are easily obtained from the diffractive cross-sections dal^ p ^ Xp /d/3, integrated over the momentum transfer t. 
In practice, one does not actually carry out the t integration of ([5]), but one rather uses the fact that the diffractive 
cross-section decreases exponentially with \t\ like e Bot . One writes: 



R , _ Q 2 /3 dal' p ^ Xp dal* p ^ Xp _ f° da[ p ^ Xp 1 dof p ^ Xp 



47r 2 a em d(3 d/3 J t . dfidt B D d(3dt 



(14) 



t=o 



with e Botmin -C 1 (in practice, t m ; m = — 1 GeV 2 ). When computing |(HJ) for t = 0, the two impact parameter integrations 
yield the factor crg/4. As already discussed, a consistent treatment of the impact parameter dependence of the dipole 
scattering amplitude T q q implies <7 2 /(4.Bd) =ttoq, which we shall use in what follows. 



Note that one could also study and F®' 4 directly. However, there is less data for those t— dependent structure 
functions, and they would not further test our model, which has the exponential decrease e Bot built in. This type of 
measurement would rather be interesting to test saturation models which feature a t— dependent saturation scale, as 
predicted in [23[ from the full BK equation. 

Let us come back to the diffractive structure functions F®' 3 . From formulae ([HI, (JEJ) and (|10p . one obtains the 
contributions from the qq final state. Using the transverse overlap function ([3|), one gets 

x P Ff((3,xr,Q 2 ) = E e // ldz ®{k))z{1-z) [{z 2 + (l^z) 2 )(z(l-z)Q 2 + m 2 )I 2 (K f ,e f ,Q s ) 

+m 2 I 2 ( Kfl e fl Q s )] (15) 

for the qq contribution to the transverse diffractive structure function. With the longitudinal overlap function (Jl|), 
one gets the qq contribution to the longitudinal diffractive structure function: 

x P Ff (P,x P ,Q 2 ) = |g^^4^ 1 dz e(4 )4 Q 2 z 3 (l-z) 3 / 2 ( K/ , e/ ,Q s ) . (16) 

In (fT^]) and (fT^|) . the functions I\ are given by 

h(K,e,Q s )= rdr.J x {nr)K x (er)N{rQ Sl x v ) (17) 
Jo 

in terms of the dipole scattering amplitude N(rQ s , xp) and of the Bessel functions J\ and K\. 



IV. THE qqg CONTRIBUTION TO THE DIFFRACTIVE FINAL STATE 

As pictured in Fig. 2, formula ([5]) is the contribution of the qq final state to the diffractive cross-section. We have 
neglected possible final states containing gluons, and in general it is justified because these are suppressed by extra 
powers of a s . However, there are two kinematical regimes for which this is not the case: the large— Q 2 limit and the 
small—/? limit. In those situations, gluon emissions are accomponied by large logarithms ln(Q 2 ) or ln(l//3) which 
compensate the factors of a s , and multiple gluons emissions should be resummed to complete formula ((5]). 

In practice, including the qqg final state is enough to describe the HERA data, and this can be done within the 
dipole picture in both limits, at leading \n(Q 2 ) [1,[2J] or leading ln(l//3) accuracy [IE, US H3, Hlj 0| i as we recall in 
this section. Note that, at leading ln(l//3) accuracy, all multiple soft gluon emissions can also accounted for in the 
dipole picture [2(1 H(| , but we shall restric this phenomenological study to the qqg contribution. 

The most popular approach is to consider the large— Q 2 limit to implement the qqg contribution [§, H3, l3l| . even 
though the experimental measurements do not reach very high values of Q 2 . In fact, the contribution of the qqg final 
state is important only for small values of /3 which, due to the finite energy available, correspond to rather small 
values of Q 2 . This is not satisfactory. In this paper, the qqg contribution to the diffractive final state is modeled in 
such a way that both the large— Q 2 and small—/? limits are implemented. 



A. The large-Q 2 limit 

At large Q 2 , the contribution of the X = qqg final state in diffractive j*p —f Xp scattering was computed in [13, H2]. 
In momentum space, the collinear gluon has a transverse momentum much smaller than Q 2 . In coordinate space, the 
scattering involves a gluonic gg dipole (see Fig. 3a): the transverse distance between the quark and the antiquark is 
much smaller than the transverse distance between the quark and the gluon. The qq pair on one side and the gluon 
on the other side form an effective gluonic color dipole which undergoes the hadronic interaction [H[ . We shall denote 
the corresponding scattering amplitude off the proton T gg (r, b; x) for a dipole of size r at impact parameter b. With 
our model for impact parameter dependence, we write 

T 99 (r, b; x P ) = 5(b) N(\r\Q s (x), x F ) (18) 



where N is the equivalent of /V but for a gg dipole. 




FIG. 3: The contribution of the X = qqg final state in diffractive j*p —* Xp scattering. Left diagram: at large Q 2 ; the quark- 
antiquark transverse distance is much smaller than the quark-gluon transverse distance and an effective gg dipole scatters off the 
proton. Right diagram: at small /3; the quark-antiquark-gluon triplet scatters after the gluon emission and the quark-antiquark 
pair scatters before the gluon emission, with a relative minus sign. In both cases only the amplitude is shown, it has to be 
squared to obtain the cross-section. 



At leading ln(Q 2 ), the qqg final state contributes only to the transverse diffractive structure function and one has 



a a s C F N c p x - 2 r 



dz 



1-2 

z 



with 



I g (a,b,c) = / rdrJ 2 (ar)K 2 (br)N(cr,xr) ■ 
Jo 



(19) 
(20) 



The computation of [24j | is a leading-twist two-gluon exchange calculation in which the gg dipole is given by 
N(rQ s ,x) — N c N(rQ s ,x)/CF in terms of the qq dipole. However this is not consistent with the use of a satu- 
ration model. For instance, N should saturate at 1, not at N c /Cf- This implies that, when using (fT9")Hj2T))) with 
N = N c N/Cf and a saturation model for N, the analysis in the literature overestimate the qqg contribution. In 
practice, this is usually compensated by using an unphysically small value for a s . 

The parametrization we shall use in this paper is N = 2N — N 2 . This relation implies the large— N c limit, and 
therefore goes well with our model for the qq dipole scattering amplitude N : it is consistent with the BK evolution 
implemented in (|13| . Numerically, this reduces the qqg contribution (with respect to using N — N c N/Cp), especially 
because the saturation scale is quite large, and therefore N is not always small. 

Finally, when computing the heavy quark contributions ccg and bbg, we replace the variable in l|19p by (3(1 + 
4m 2 /Q 2 ). This substitution, which modifies only the small— Q 2 results, is necessary in order to insure that there is 
no qqg contributions when the final state is such that Mx = 2m/ (in practice, such a substitution doesn't make a 
difference for the light quarks). 



B. The small-/3 limit 



At small 0, the contribution of the X = qqg final state in diffractive "f*p — > Xp scattering was computed in many 
studies [IH, IH, H3, HI, Hi!- In coordinate space, denoting x the transverse position of the quark, y that of the 
antiquark and z that of the gluon, the diffractive scattering is expressed in terms of (see Fig. 3b) 



[T, 5fl (x, y, z; a*) - T, 9 (x, y; x P )] 2 = S 2 (b) \n^(\v'\Q Si |r - r'\Q s , x r ) - N(\r\Q s , x r ) 



(21) 



In the left-hand side, the virtual contribution T q q represents the scattering of the quark-antiquark pair, before the 
gluon emission. Within our model for the impact parameter b=(x+y)/2, one has T q q(x, y; xp) = S(h)N(\r\Q Sl Xf) 
where the dipole size is naturally r = x— y. The real contribution T q q g represents the scattering of the quark-antiquark- 
gluon triplet, after the gluon emission. In the right-hand side, we factorized the impact parameter profile and wrote 
T qq - g (x,y,z;x P ) = S(b)N( 2 1{\r'\Q s ,\r - r'\Q s ,x r ) with r' = x-z (and r - r' = z-y). 

In the context of the BK evolution implemented in (fT3|) . the link between ./V*- 2 ) and N comes from the fact that the 
scattering of the qqg triplet is equivalent to the scattering of two dipoles with sizes r' and r— r' (a dipole emitting a 
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FIG. 4: The contribution of the qqg final state to the transverse diffractive structure function at (3 = as a function of Q 2 . 
The full lines show the exact result F^? a \ LL (^j^ while the dashed lines show the leading ln(Q ) result ^t' s |lz,(q 2 )- Different 
sets of curves are for different values of = 0.01, 0.001, 0.0001, 0.00001, from bottom to top. As Q 2 increases the two results 
get closer, but they coincide for only very large values of Q 2 . 



soft gluon is equivalent to a dipole splitting into two dipoles). Therefore our model for is 

N^(\r'\Q s , \r-r'\Q s ,x r ) = N(\r'\Q s ,x F ) + N{\r-v'\Q s ,x P ) - N(\r'\Q s ,x F )N(\r-r'\Q s , x P ) . (22) 
At leading ln(l//3), the contribution of the qqg final state to the transverse diffractive structure function is 

xr F^\ LL(im (x P ,Q 2 ) = ° Fa f 2ao r rdr [ dz $ T (z, r; Q 2 )A(r, x r ) (23) 

with 

A(|r|,a*) = J d 2 r' r , 2 ^ 2 _ r/)2 [N(\r'\Q s , x r ) + N(\r-r'\Q s , x P ) - N(\r\Q s , x v ) - N(\r'\Q s , x r )N(\r-r'\Q s , x P )f . 

(24) 

It is independent of j3 because the structure function F D ' 3 ~ (3 da 1 p ^ Xv /df3 picks up the coefficient of ln(l//3) in 
cr 7 p~* x p. Also, the overlap function is <&t because in the leading ln(l//3) approximation, the final state mass Mx is 
fixed only by the soft gluon longitudinal momentum, and therefore transverse sizes are the same in the amplitude and 
the complex conjugate amplitude. 

Note that in (|23| . the impact parameter integration J d 2 b <S* 2 (b) yielded a factor cr°/4. If the b— profile was a theta 
function as assumed in [28|], the qqg contribution would be a factor of 2 higher. In what follows, to numerically 
compute A(r, xp), we use the clever change of variables introduced in (2e| that we recall in Appendix A. 



C. The model for x P F™ 3 

The usual approach to implement the qqg contribution is to use formula (fTT)|) . but as we shall see, this is not correct 
for small values of (3. Let us consider the qqg contribution for (3 = 0. By definition, the correct result is F^ q9 \i n ^i^ 
given in formulae ([2l?lH[2"4")) . By contrast, the small—/? limit of the leading ln(Q 2 ) contribution F^ q9 \ LL (Q2j is 

f°° dr 2 
/ — J 2 (kr)(2N(rQ s ,x P )-N 2 (rQ s ,x P )) 
Jo r 

(25) 

where we have used (fT9|) with K 2 {x) = 2/x 2 for x — > 0. Formula ([25]) shows that, after rising as (3 decreases, the 
diffractive structure function goes to a constant. This constant is different from the correct result (f2"3")Hfl?4")) . except 
for very large values of Q 2 , for which Fj,"\ll(Q 2 ) is correct by definition. And indeed, if Q 2 ^>Q 2 , the two formulae 



XV r T \LL(Q 



2)(/3 = 0,x P ,Q 2 ) 



C F N c a s a 
12tt 4 



dk 2 



In 



k 2 



coincide to give 



x P F^ 9 (P = 0,x P ,Q 2 »Q 2 ) = C I N0la 1 ^ e , ln ^ J™ J [ 2N {f , x r ) - N 2 {f , x r )} 2 . (26) 

This is shown analytically in Appendix B. In Fig. 4, we compare formulae ([2"3 |H[2"4"|) and (|25|) as a function of Q 2 and 
for different value of Xp. One sees that when Q 2 increases, the ratio between the two results gets closer to one, but 
that limit is only reached for very large values of Q 2 not shown in the figure. For the values of Q 2 in the HERA range, 
the actual result is smaller than the leading ln(Q 2 ) one by a factor of about 0.6. 

In order to have the correct qqg contribution for small values of (3, we shall use the following model: 

Xr F^{(5,x^Q 2 ) = xpF T 99 \ LL{Q2) (f3,x r ,Q 2 ) JP LL{ Z ){ 7 Q \ W 

obtained from formulae (PjH|2"0|) . pJjHjMj) and (J2SJ. It is such that Ff 9 = F^ 9 \ LL ^ at large Q 2 and Ff 9 = 

Ft"\ll(i/i3) a t small f}. In the small— Q 2 and large—/? region, the qqg contribution may not be correctly described. 
However in this case, the diffractive structure function is dominated by the qq component, and the qqg contribution 
is not relevant. 

Finally we point out that our implementation of the qqg contribution is parameter free, the only uncertainty being 
related to the value of a s . The average value of Q 2 in diffractive measurements at HERA is about 10 GeV 2 , therefore 
we choose a s = 0.25, which corresponds to such a scale. 



V. DESCRIPTION OF THE HERA DATA 



The HI and ZEUS experiments at HERA have measured the diffractive cross section for the process ep — > eXp, 
tagging the proton in the final state. After integrating the squared momentum transfer dependence from t m i n = 
— 1 GeV 2 to i = 0, the data are presented in terms of the reduced cross section a®' 3 ((3,xp,Q 2 ) : 

dxVWW'^^K V+ ^J r (/? ' XP ' 0) ' ~ + 2-2y + y 2 ^ ' (28) 

with y = Q 2 /(sx) where -^5 = 318 GeV is the total energy in the e—p collision. We shall call the corresponding data 
sets the LPS [H| (ZEUS) and FPS ^ (HI) data. 

The HI and ZEUS experiments have also measured the diffractive cross section for the process ep^eXY, selecting 
events with a large rapidity gap between the systems X and Y in case of HI (33|, and using the so-called Mx- 
method in case of ZEUS [34| . Y represents the scattered proton, either intact of in a low- mass excited state, with 
My < 1-6 GeV (HI) or My < 2.3 GeV (ZEUS). The cut on the squared momentum transfer t at the proton vertex, is 
again t > — 1 GeV 2 for both experiments. We shall call the corresponding data sets the FPC [HI (ZEUS) and LRG 
[H (HI) data. 

Because they include events in which the proton has broken up, the cross-sections measured for the process ep^eXY 
are larger than the one measured for the process ep— >eXp. Also, because HI and ZEUS measurements are performed 
with different My cuts, the ZEUS cross-section is bigger than the HI cross-section, for which the proton-dissociative 
events are more reduced. However, within the kinematical ranges of the measurements, it seems that the differences 
are constant factors: the FPC and LRG data points can be converted to the FPS-LPS ones by dividing the cross- 
sections by 1.45 and 1.23 respectively [33 ! , 34]. Note that it is the FPS-LPS data that correspond to our definition of 
diffractive events and to our formulae, as the proton should truelly escape the collision intact. 

In our model, the reduced cross section a® ,3 ((3,xp,Q 2 ) is given in terms of the diffractive structure functions by 

x v a^ = xpF 99 + xpFf 9 + - 2 ~ 2y xpFf . (29) 

2~2y + y z 

Using formula (fl"5|) . (fT!)|) . (|27[) . and the dipole scattering amplitude (jX3J> , we obtain a parameter free calculation for 
xpa®' 3 that we can compare to the data. Diffractive DIS measurements are sensitive to the saturation regime of 
QCD only for small values of xp, therefore we shall only consider experimental data which feature xp < 10~ 2 in our 
comparisons. Note that we do not include any qqg contribution to the longitudinal structure function: for small values 

D 3 

of f3 it could be sizeable, but for kinematical reasons small (3 is associated with y close to 1, in which case F L ' does 
not contribute to tri 3 ' 3 . 
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FIG. 5: x r <T?' 3 (/3,xi.,Q 2 ) as a function of x r for different values of (3 and Q 2 . The HI (LRG) and ZEUS (FPC) diffractive 
data are compared to the predictions of our model, and the xr range is restricted to a* < 0.01. In this figure, the HI data are 
unchanged and it is our predictions which are multiplied by the factor 1.23. The ZEUS data are multiplied by 0.85 in order 
to convert them to the HI My range and the bins centers have been shifted to the HI values using a parametrisation given 
in [34)]. Only the statistical part of the uncertainty is shown for the data points on this plot. The shape of the curves in the 
(3 = 0.9 bins is due to the fact that the contribution of F L ' to the reduced cross-section is important for large values of (3. 




To estimate the quality of our description, we performed the following \ 2 computations, adding statistical and 
systematic uncertainties in quadrature. Within the LPS+FPS (ep^ eXp) data sets, 76 points pass the ccp<0.01 cut 
and we obtain x 2 /points = 0.80. When comparing to the 4 data sets, with the proper renormalizations for the FPC 
and LRG (ep— > eXY) measurements, 343 points pass the xr < 0.01 cut and we obtain x 2 /points = 1.28. This is a 
quite good description, considering our predictions are parameter free. As an illustration, Fig. 5 displays a comparison 
of our predictions with the FPC+LRG data. Wc also checked the agreement with the charm contribution to a®' 3 
using the few points available (35[, one obtains % 2 /points = 0.68. 

Note that, if the running coupling a s (Q 2 ) is used in the qqg contribution (instead of imposing a s = 0.25), the 
description of a®' 3 is also very good with similar values of x 2 /points : 0.63 when comparing to the LPS and FPS 
data, and 1.20 when comparing to the four data sets. In this case, the additional cut Q 2 > 1 GeV 2 is used in order to 
keep the coupling reasonably small (this only removes 15 of the LPS points). We also noticed that it is possible to 
obtain a description of equal quality without the correction F. '^p 9 / p) I F!^ 9 \ll^q2) in the qqg contribution (|27|) . if 
one imposes an unphysically small value for the coupling: a s = 0.15. 

Finally, Fig. 6 shows predictions for the longitudinal diffractive structure function F®' 3 ((3,xp,Q 2 ), which in our 
approach in obtained from the qq contribution F L q (see formula (|16p ). 



VI. CONCLUSIONS 



We presented a new description of HERA diffractive deep inelastic scattering data. It uses the parametrization for 
the dipole scattering amplitude obtained in |IH ■ This is an extention of the IIM saturation model which contains 
heavy-quark contributions. Contrary to previous studies, it features a saturation scale Q s (x) that stays above 1 GeV 
for :r = 10~ 5 , rather than dropping by a factor of 2 after including heavy-quark contributions. Our description of the 
data is parameter-free and, in the regime xp<0.01 in which saturation is expected to be relevant, it features values 
of x 2 /points of order 1. 

Let us recall the improvements that our model brings with respect to previous approaches. 

• Instead of considering the dipole cross-section ao and the diffractive slope Bo as independent quantities, we 
use the relation pip . It results from a consistent treament of the impact parameter dependence of the dipole 
scattering amplitude. 

• In the large— Q 2 limit, the qqg contribution to the diffractive final state is described by a gg dipole. It is related 
in terms the qq dipole in such a way that it is consistent with the BK evolution implemented in (|13p . 

• The qqg contribution to the diffractive final state is modeled in such a way that both the large— Q 2 and small—/? 
limits are implemented (|27|) . This allows to have a good descrition of the data with a consistent value of 
a s = 0.25. 

F) 3 

• Our predictions include correclty the contribution of the longitudinal diffractive structure function F L ' : we 
predict a?> 3 , not F?> 3 = F°' 3 +F° < 3 . 



• When comparing our predictions to the experimental data, the comparison is made with the ep — > eXp data 
unchanged and the ep—>eXY data renormalized, and not the opposite. Our definition of diffractive events (and 
our formulae) are such that the proton truelly escape the collision intact. 

Of all the possible descriptions of diffractive DIS data (for a global analysis, see [Ifl), the dipole picture is the one 
which is adapted to study the physics of parton saturation. Having a consistent saturation model to describe hard 
diffraction in e — p scattering represent a good foundation for further works. On the phenomelogical side, studying 
hard diffraction in e— A becomes of interest [37j . as it is an ideal process to investigate the saturation regime of QCD 
that could be explored at a future electron-ion collider. On the theoretical side, many new developements improved 
our understanding of the QCD non-linear evolution, and their consequences are to be investigated [20I |38| . 
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APPENDIX A: NUMERICAL COMPUTATION OF A(\r\,x r ) 

In this Appendix, following [28^ . we show how to transform the expression (see formula (|24jl) 

A(\r\,x v ) = / d 2 r' * [N(\r'\Q s ,x r ) + N{\r-r'\Q B ,x r ) - N(\r\Q s ,x r ) ~ N(\r'\Q s , x r )N{\r-r'\Q s , x P )} 2 
J r'-^r — r' y 

( A1 ) 

in order to estimate it numerically. Writing the two-dimensional integration in the complex plane and introducing 
S = 1 — N, we obtain (with r= |r|) 



dzdz 



A (r,Xf) = I — . 2 [S(\z\rQ s ,x P )S(\l-z\rQ s ,x F ) - S(rQ s ,x P )] 2 . (A2) 
J 2\z\ z \l-z\ z 

We then follow the following procedure. 

• \z\ and |1 — z\ are invariant by symmetry with respect to the real axis so one can multiply the integral by 2 and 
restrict ourselves to the upper part of the complex plane. 

• For \z\ < 1, let us change the variables into u = \z\ G [0, 1] and v = ^—^U^—^ G [0, 1], this implies 

dzd - = 8\z[j^-z\dudv = 4:u(l-u + 2uv)y/v(l-v)(l+uv)(l-u + uv) . (A3) 

\z—z\ 

• For \z\ > 1, let us change the variables into u = l/\z\ G [0, 1] and v — 1 ^ z ^ 2 IP' 1]j this implies 

, , 8|z| 3 |l — zldudv \l — z\, , , , 1— —, r-, r ,, , 

dzdz^-^ ^ J—— -\z-z\ = A{l/u-l + 2v)yjv{\-v){\+uv){\-u + uv) . (A4) 

\z z\ \z\ 



One gets 

, rl f 1 2dudv 

A{r, x r ) 



Iq Jo u(l — u + 2uv)y/v(l — v)(l + uv)(l — u+uv) 

{u 2 [S(rQ s /u,x r )S((l/u-l + 2v)rQ s ,x P ) - S(rQ Sl xr)] 2 

+ [S(urQ s ,xp)S((l-u+2uv)rQ s ,x v ) - S(rQ s ,x v )] 2 } (A5) 

which is easy to evaluate numerically, as this features only integrable singularities. 



APPENDIX B: THE SMALL—/? AND LARGE-Q 2 LIMIT OF x v F* 



In this Appendix, we show that the small—/? limit of the leading ln(Q 2 ) result F!j? a \ LL iQi^ and the large— Q 2 limit 
of the leading ln(l//3) result F^ q& / (i) coincide. The first case has been derived in the text and we obtained formula 



We now show how, for Q 2 ^>Q 2 , formulae ([2"5)Hf2"4")) give the same result. 
The starting point is the transverse overlap function 



®T(z, r] Q 2 ) 



2tt 2 



^e 2 ([z 2 + (l-zf]s)K 2 {re f ) + m 2 K 2 (re f )) , 
/ 



(Bl) 



where one can neglect quark masses with respect to Q 2 . The z integration of (|B1[) can be done in the two limits 
rQ<Cl and rQ^>l : using the Mellin representation of K 2 (x), one gets 



dz[z 2 + (l-z) 2 }z(l-z)K 2 pyjz(l-z) 



c+ioo 



d7 _ 27 r(7-i)r( 7 )r( 7 +i)r(2- 7 )r(4- 7 ) 



2m' 



r(6-2 7 )r( 7 +i/2) 



l<Re(c)<2 



2 ( 4/p 4 for p > 1 
1/p 2 for p < 1 



(B2) 



To obtain the second equality, we used the fact that in the pS>l case (p<Cl case), the dominant contribution to the 
7 integration comes from the single pole at 7 = 2 (at 7=1). One can then write 



nir i dz* T (z,r;Q 2 )f(r) = ^^Y. e2 f 

71 f 



2 ' Q dr 



f{r) 



Q 2 



dr 



f{r) 



2/Q 



(B3) 



Note that in [28j |. a similar estimation was obtained. By contrast, the replacement K\{p)—>Q(l — p)/p was used and 
as a result, the normalization factors were not under control. 

To complete the calculation, we need to compute f(r) = A(r, xp), in the two limits rQ s <Cl and rQ s 3>l (see [28|): 



A(r,xp) = 2nr 2 



J [2JV(f, x P ) - N 2 (f, xp)] 2 for rQ s « 1 



A(r,x r ) = 27rln(r 2 Q 2 )[l - N{rQ s , x P )} 2 for rQ s »l . 



(B4) 
(B5) 



The first line is obtained from configurations in which the qq pair is small and is well separated from the gluon. 
When using those results in (|B3|1 by dividing the r integration region in three domains, one sees that the dominant 
contribution comes from the region r e [2/Q, 1/QJ. It is enhanced by the collinear factor ln(Q 2 /Q 2 ) : 



xp F« 9 (/3 = 0,Q 2 »Q 2 ) 



C F N c a s Q 2 <T 
6tt 4 



E4M|) /„ ^[2N(f,xp)-N 2 (f,x v )Y 



dr 



(B6) 



This formula is identical to formula ([2^)1 . 
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